function [TUteJTPAlb,TUteJTPAub] = JTPA_UNION_TU(n0,n1)

load JTPA_CATE_UNION_INIT;

% Precompute crossproduct matrices
MZZ = Z'*Z;  
MZY = Z'*y;  

% Targeted undersmoothing
nrep = n1-n0+1;
TUteJTPAlb = zeros(nrep,10);
TUteJTPAub = zeros(nrep,10);
for ii = n0:n1
    fprintf('Iteration: %d \n',ii);
    [TUteJTPAlb(ii-n0+1,:),TUteJTPAub(ii-n0+1,:)] = fsel_linear_big(Z,y,XteUnique(ii,:)',usepl,MZZ,MZY,10);    
end

filename = strcat('JTPA_CATE_UNION_TU_',num2str(n0),'_',num2str(n1));
save(filename,'TUteJTPAlb','TUteJTPAub','n0','n1');

